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Abstract. In this paper we present a new family of solutions to the singularly perturbed 
Allen-Cahn equation a^Au + «(! — u^) = 0, in C R'^ where N = 3, Q, is a smooth 
bounded domain and a > 0 is a small parameter. We provide asymptotic behavior which 
shows that, as a —>■ 0, the level sets of the solutions collapse onto a bounded portion of a 
complete embedded minimal surface with finite total curvature that intersects orthogonally 
dfl of the domain and that is non-degenerate respect to fl. We provide explicit examples of 
surfaces to which our result applies. 


1. Introduction 

1.1. Preliminary discussion. In this paper we study the singulary perturbed boundary 
value problem 

du 

a'^Au + u{l — u‘^) = 0, in — = 0, on dD. (1.1) 

on 

where a > 0 is a small parameter, H C IR'^ is a smooth bounded domain and n is the inward 
unit normal vector to dfl. 

Solutions to dni) correspond exactly to critical points of the Allen-Cahn energy 

Jc{u) := f + ^(1 - S 

Equation (11.111 arises for instance in the gradient theory of phase transition when modelling 
the phase of a material placed in D or when studying stationary solutions for bistable reaction 
kinetics, see [2]. 

Observe that r = ±1 are global minimizers of Ja, representing stable phases of two different 
materials isolated in H. 

We are interested in solutions in solutions u connecting the stable phases ±1. As described 
in [28], solutions of this type are expected to have narrow transtition layer from —1 to -|-1 
with a nodal set that is asymptotically locally stationary for the perimeter functional. To be 
more precise in |28j . the author showed that a family of local minimizers {uo,}a of Ja with 
uniformly bounded energy must converge in L^(H), up to a subsequence, to a function r*, 
where 

u* =XA- Xn-A 

where xe is the characteristic function of a set E and A C D minimizes perimeter. In this 
case, as a —)• 0 

Ja{ua) -A Per{A) ^\w'f + ^(1 - (1.2) 
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where w is the solution of 

w" + w{l — w^) = Q, in R, in'> 0, r(;(±oo) = ±1 (1-3) 

which is given by the explicit function 

t 

Roughly speaking, the above assertion means that the level sets of Ua converge to 5A as 
a —7- 0. This result provided the intuition that ultimately lead to important developments in 
the theory of T—convergence and put into light a deep connection between the Allen-Cahn 
equation and the theory of minimal surfaces. We refer reader to [i El IH ED [33] for related 
results and stronger notions of convergence. 

The connection between the Allen-Cahn equation and the theory of minimal surfaces has 
been explored in order to produce nontrivial solutions of (11.11) , but the general understanding 
of solutions to these equation is far from being complete. In this regard it is natural to ask for 
existence and asymptotic behavior of solutions to o in the general setting. For the case of 
minimizers we refer the reader to [DllolEaES] and references there in. We also remark that 
local minimizers in convex domains are the constants ±1, see [271|6|. 

In low dimensions, for instance, N = 2 Kohn and Sternberg in [T2|, using a measure 
theoretical approach and taking advantage of the aforementioned intuition, constructed local 
minimizers Ua to (|l.lll with interfaces collapsing onto a fixed minimizing segment Tq inside 
that cuts d^l perpendicularly and satishes (11.21) . 

In 1 13] a situation similar to that one described in [T2| is consider but lifting the minimizing 
assumption on Tq, for nondegeneracy of this segment respect to the domain. Nodegeneracy 
of To in this case is stated as 

/Fo + All - iToli^oi^i ^ 0 

where Kq,Ki are the curvatures of cAI at the points where Tq cuts cAI orthogonally. 

This geometrical condition is equivalent to the fact that the eigenvalue problem 

h" = Xh, in(0,/), Koh{0)+ h'{0) =0, Kih{l) - h'{1) = 0, Z = ITqI 

does not not have A = 0 as an eigenvalue. The author also provides a full description of the 
solutions which can have Morse Index either one or two depending on the sign of Kq and Ki. 

Later this construction was generalized in under the same geometrical setting described 
in m, but for multiple transitions that in the limit collapse onto the segment Tq. There the 
transition layers at main order interact exponentially respect to their mutual distances giving 
rise to Toda system of odes. 

In dimension N = 3, Sakamoto in [35| constructed solutions to (El) having a narrow 
transition through a planar disk that cuts orthogonally the boundary of the domain and 
is nondegenerate in a suitable sense. The author also provides a characterization for this 
nondegeneracy in terms of spectrum of the dirichlet to neumann map of the planar disk. As for 
higher dimension in the setting of manifolds, Pacard and Ritore in |30j constructed solutions 
having a narrow transition along a codimension one nondegenerate minimal submanifold. 

In the the spirit of the results mentioned above, we also want to refer the reader to mm 
ESKISI dealing with similar results for the inhomogeneous allen-cahn equation and [371 (381 [39] 
for a semilinear elliptic problems where resonance phenomena appear. 


w{t) = tanh 


t e R. 
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The underlying geometric problem is the existence of minimal surfaces inside the domain 
n intersecting the boundary dQ orthogonally. This problem, in a general three dimensional 
compact Riemannian manifold, has been completely settled in a recent paper by M. Li [26]. 
For earlier results in this direction we refer to Fraser-Li [19] and Fraser-Schoen |17j . The 
uniqueness of the critical catenoid in a ball is proved by Fraser-Schoen m- 

1.2. Main result. Our goal in this paper is to generalize the results in [131135j respectively 
by taking N = 3, and a more general class of minimal surfaces for limiting nodal set. 

In order to state our main result, let M be a complete embedded minimal surface of 
finite total curvature in For over a century there were only two known examples of 
such surfaces, namely the plane and the catenoid. In [7] [ 8 ] Costa gave the first nontrivial 
example of such surface with genus one, being properly embedded and having two catenoidal 
connected components outside a large ball sharing an axis of symmetry and another planar 
component perpendicular to this axis. Later this construction was generalized in [SUES] to 
surfaces having the same look as the Costa’s surface far away but with arbitrary genus. We 
refer the interested reader to [25[ [32] and references there in, for related results and further 
generalizations. 

It is known that M is orientable and — M has exactly two connected components namely 
S~^ and S~, see [23]. We set by u : M —)• 5^ the fixed choice of the unit normal vector of M 
pointing towards S~^. 

For X = {xi,X 2 ,X 3 ) G we denote 

r = r(x) = ^x‘l + X 2 - 

It is also known that for some large but fixed Rq > 0, outside the infinite cylinder r > Rq, 
M decomposes into finite connected components, say Mi ,..., M^, which from now on we will 
refer to as the ends of M. For every k = 1,... ,k, there exist a smooth function = F^iy') 
such that 

Mk = {(?/',ya) G : r{y',yi) > Rq, 2/3 = Fk{v')] 
where Fk has the asymptotic expansion 

Fk{y') = ak'^og{r) + hk + bik^+0{r~^), as roo (1.4) 

for some constants ak,bk,bik satisfying 

m 

(ll ^ 0,2 — ' ' ' — Cik, ^ ^ Rfc — 0 

k=l 

and relation (jl.4|) can be differentiated. In this case we set the connected component of 
R^ — M to be the one containing the axis X 3 which corresponds to the axis of symmetry of 
the ends Mi,..., M^. 

We denote by u : M —)• the unit normal vector to M pointing towards S'"*" and we take 

Fermi coordinates near M 

x = y + zv{y), yeM, \z\ < rj + 6log{2 + r{y)) 

for some 77 ,5 > 0 small. From now on 6 will represent a small, but fixed number, independent 
of q; > 0. Observe that also that R^ is spanned by the moving frame Ty { M ) 0 Ru { y ) and that 
z corresponds to the signed distance to M, i.e 

\z\ = dist{x, M), X = y + zv{y) 
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for every y G M and z small. 


Next, consider a smooth bounded domain such that 

(i) . n contains a portion of M which we denote by M. 

(ii) . n — M has two connected components which, abusing the notation, we naturally 

denote as S’"*", S~ with the same convention as before. 

(iii) . dVL n M = n dVL where for every fe = 1,..., m, Cfc := Mk n 90 is a smooth 

closed simple curve. 

Notice that the 90n7W consists of nonintersectiong closed curves, since the curve Ck C Mk- 

Following [I2l[l3l[9] in order to produce solutions of (II.Ih we need some kind of criticality 
and nondegeneracy assumptions on M. respect to O. To do so, let us introduce A_m the 
laplace-beltrami operator of M, and \A_m\ the norm of the second fundamental form of M 
and let us consider now the following eigenvalue problem 

dh 

Aj^h + \AM\‘^h = \h, in A^, — — \- K{y)h = 0, on (1.5) 

where r represents the inward unit normal direction to A4 and K{y) is given by 


K{y) := 


dn \ 


y £ M. 


where we recall that n is the unit normal vector to dVl. We remark that since A4 is a minimal 
surface it follows that 

with Kj^ being the gaussian curvature of A4. 


Our crucial assumptions on A4 are the following. 

(I). A4 cuts orthogonally dQ along the curves (7^, for /c = 1,..., m. 

(II). The eigenvalue problem (|1.5h in does not have A = 0 as an eigenvalue. 


As stated in [39], from assumption (I) it follows that r and n must be parallel along every 
curve Ck and consequently these curves must be geodesics in 911 in the direction of ly since 
their the normal vectors in 9fl are parallel to n. Therefore, the quantity K{y) corresponds to 
the geodesic curvature of 911 in the direction v{y) for y £ Ck- 

Our main result is the following. 


Theorem 1. Assume conditions (i)-(ii)-(iii) and (I)-(II)- Then for every a > 0 small enough 
there exists a solution Ua to dEU) such that 


Ua{x) = w 

for X = y + zuiy) £ Af = {dist{-, M) < rj} nil, where the function h solves at main order the 
boundary value problem 


dh 

+ |A_a 4 |^/i = 0, m A4, — \-n{y)h = 0{a)^ on 9A4 


While outside M, 


Ua(x) = 


dr 

+1, X e 5+ 

— 1, X £ S~ 


( 1 . 6 ) 


(1.7) 


as a —0. 
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Theorem [T] provides a rather general description of an explicit family of minimal surfaces 
for which the construction applies and allows us to consider more involved examples of the 
one presented in |35j . We also remark that if our original surface and the domain have axial 
symmetry one can reduce our developments to this setting and condition (II) can be recasted 
as requiring that A = 0 is not an eigenvalue, for problem (11.51) in 

The paper is organized as follows. In section [2] we present briefly invertibility theory of the 
operator descibed in (II. 6 p with robin boundary conditions and we also disscuss some example 
where our result apllies. Next in section [3l we present the geometric framework we will use 
to set up the proof of Theorem [TJ In section 0] we construct an accurate approximation of the 
solution to our problem and then in section [5] we sketch the proof of our main result. The 
final section is devoted to present detailed proofs of lemmas and propositions used in section 

El 

Acknowledgments: The research of the first author was supported by the Grant 13- 
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2. JACOBI OPERATOR WITH ROBIN BOUNDARY CONDITIONS 

In this part we consider the equation 

dh 

AMh+\AM\‘^h = f, — + I{y)h = 0, on dM 

Using Fourier decomposition and the nondegeneracy assumption (II), it is straight-forward 
to verify that for any / S L^(AI), there exists a unique solution h S satisfying 

\\h\\w2,2(^M) ^ C'II/IIlRAJ)- 

If p > 2 and / € LP{M.) then by standard regularity theory h G n 

with a priori estimate 

||h||* = \\D'^h\\Lp(^M) + + \\h\\L°°{M) < C\\f\\LP{M)- 

Directly from this it follows that for any / G LP{Ai)^ g G LP{dM) 

Bh 

AMh + \AM\'^h = f, in Ad, —+I{y)h = g, on dM 

or 

has a unique solution satisfying 

||/i||* < C (||/||lp(aj) + WaWiPidM)) ■ 

2.1. Examples. In this part we discuss some particular situations where our theorem applies. 
Let M be the catenoid in IR^ parameterized by the mapping 

■= (\/l + y2 cos 6l,Vl-ky2 sine, log (^y-R Vl-ky2^^ , y G R, 6lG(0,27r) 

which provides coordinates on M in terms of the signed arch-length of the profie curve and 
the rotation around the xa-axis, which in our setting corresponds to the axis of symmetry of 

M. 

The unit normal vector to M pointing towards S'"*', is given by 

p(y,^)= „ (-cos g,-sing, y), y G IR, 6» G (0, 27r) 

VI+F 
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and we consider the Fermi coordinates 

^(y, d,z)=Y (y, e) + zv{y, 9) 

which define a change of variables for instance on the neighborhood of M, 

^ = ^{y,^) + ziz{y,9) : |z| < r/+iln(l+y^)| 
for some fixed and small r] > 0. 

Now we assume that fl is axially symmetric. Since M has two ends and due to the axial 
symmetry, dQCiM = C 1 UC 2 where Ci,C 2 are parallel, nonintersecting circles parameterized 
respectively by 

Ue):=Y{y„e), 0G(O,27r), z = l,2 

for some fixed y^ < y 2 - 

To describe dit close to the circles Yi(6), we can assume the existence of two smooth 
functions 

Gi,G 2 : ^ GiiO):=yi, i = 1,2 

so that the two systems of coordinates 

X^{e,z):=Y{Giiz),6) + zi^{G^{z),e), 0G(O,27r), \z\<7], i = l,2 (2.1) 

describe the set 

dit n AA = {x € dit : X = Xi{6, z), \z\ <7], 0 € (0,27r), i = 0,1}. 


In order to explain geometrically conditions (I), (II), first recall that Ai = Q Ci M and 
consider a function h G G‘^{Ai) such that ||h||( 72 (_v)) < rj. A normal deformation of Ai within 
n can be described by the coordinate system 

Yh{y,9) ■.= Y {y{y,e),e) + h{y,e)u{y{y,e),d) , yi < y < y 2 , 6» G (0, 27r) (2.2) 


where 
y(y, 6 ») := 


G2{h{y,e))-G2{h{y,e)) 


{y-yi) + Gi(h{y,e)), yi<y<y 2 , 6»G(0,27r) 


y 2 -Yi 

Denoting Aih '■= Yh{[yi,y 2 \ x (0,27r)) and det gh its respective induced metric, with the 
convention that go is the induced metric of Ai. The area functional oi Aih is computed as 

f27r i-y^ _ 

^Qh - ! I Vdet ghdyde. (2.3) 


A{A4h):= [ IdA. 

JMh -^0 Jyf) 

This area functional is of class and its first variation around Ai is given by 
'd,G2iO) - d,Gi{0) 


DA{Ai)[h] = - [ 

Jc 


C 1 UC 2 


y 2 -Yi 


where Hj^ is the mean curvature of Ai. 


(y-Yi) + dzGi (0) ) h{y, 6)dsg^ + / Hm h dAg^ 
) Jm 

(2.4) 


We notice from (12.41) that Ai is critical for the area functional (12.31) if 

Hm = 0, 5,Gi(0) = 5,G2(0) = 0. (2.5) 

Therefore condition (I) is equivalent to saying that Ai is critical for the functional (12.3p 
respect to normal perturbation of Ai. 


Since A4 is a minimal surface, condition Hm = 0 is automatically satisfied. Consequently, 
the criticality of Ai reduces to saying that Ai must cut dO. orthogonally. 
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Assuming condition (I), the second variation of the area functional around A4 is given by 
the quadratic form 

DU(M}lh,h] := [ d,,GiiO)h^{y)dsg, + [ dAg, 

Jci J M 

and stability properties of Ai respect to are analyzed through the linear robin boundary 
value problem 

dh 

+ |A_A/(p/i = A/i, in A4, - -h i9^2Gj(0)/i = 0, ondM. (2.6) 


We make the important remark that, under assumption (I), re = r* where r* is the inward 
unit tangent vector of the profile curve of the catenoid along the circle Cj. 

Also observe that Ki := dzzGi{0) corresponds to the curvature of the integral curve 

of in the direction of re along the circle Cj. 

From the axial symmetry, nondegeneracy of Ai reduces asking that the only solution to 
the boundary value problem 

dyyh + + (1 + y2)2 ^ = Q’ yi<y<y2, ^(y*) + (-i)*-?^i^(yi) = o. 

is the trivial one. 


The linear equation 

Am^ + |AAfph = 0, in M (2-^) 

has two axially symmetric entire solutions 

zi{y) = y-e3, Z 2 {y) = y ■ v{y), yeM 

corresponding respectively to the invariances of the entire catenoid M under translations 
along the vertical axis and dilations. We refer the reader to section 4 in [T] and m for full 
details. 

We directly check that in coordinates y = Y{y,6) G M 

^i(y) = -7=^^ ^ 2 (y) = log(y + Vl + y2) - 1, y G K 

Vi + y^ Vi+F 

from where we observe that zi is odd and Z 2 is even. 


We remark also that zi,Z 2 are strictly increasing in the variable y G (0,oo), since 

y log (^y + Vl + y2 

1+y^ (l+y2)l 

and also Z 2 changes sign at only one point y = yo > 0 and 2^2 (y) < 0 for —yg < y < yo- While 
zi changes sign once at y = 0. 


dyZiiy) = 


1 


(1 + y 2)2 




dyZ2{y) 


y G K 


Using basic theory of odes and following the developments from section 4 in [T], we find 
that A = 0 is not an eigenvalue of (j2.6p if and only if 


det 


dyZiiyi) + Kizi{yi) 5y 2:2 (y 1 ) + All 2:2 (y 1 ) 
dy^i (y 2 ) - K2Z1 (y2) dyZ2 (y2) - K2Z2 (y2) 


^ 0 . 


( 2 . 8 ) 


In particular, we directly check that in the case that cAI is almost flat, i.e Ki = K 2 = 0, 
Ai is nondegenerate. We also remark that condition (12.81) is clearly invariant under dilations. 
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, ^2 2 2 
0 := ^ + ^ + ^ = l 


Consider the case 


an ellipsoid of revolution, where a, 6 > 0 and A4 being and even catenoidal portion inside n. 
In this case from the disscussion above, we can describe dfl near A4 as the level set 


4 f Vi + y2 - 


+ i (iog(y + Vi + y^) + 2 


Vi + y^ 


62 




= 1 . 


Using implicit function theorem we find that 


d,G2(0) = -d,Gi(0) = 


+ y') 


y 




so that for to be critical respect to the ellipsoid II we need to choose the point y satisfying 

and from the monotonicity of 2^2 (y) we see that once the ellipsoid has been fixed there is 
exactly one catenoid that cuts the boundary of the ellipsoid perperdicularly. This is our 
candidate of critical and nondegenerate minimal surface inside the ellipsoid. 

In order to compute the geodesic curvatures we apply again the implicit function theorem 
to find that 

K 2 = Ki = -h 

a 


In the case H with R = a = b, we observe that Ai is the so called critical 

catenoid. This setting appeared in m as solution of an maximization problem for the first 
steklov eigenvalue of the dirichlet to neumann mapping in bounded domains. 

Using the same notation as above, one can verify that Ki = K 2 = ^ and the 

nondegeneracy condition translates into saying that f ^ which holds true since y ~ 1.5. 

Concerning stability issues let us consider the quadratic form 

q{h,h):=- [ k{y)h^dsg,+ [ {\Vh\^ -\AM\^h^)dAgg, h G H\M). 

JdM Jm 


We first stablish conditions Ai to be minimizer of the area functional. 


Proposition 2.1. Assume z is a smooth positive solution to the linear equation (12.7p in an 
open set of M, containing Ai. For every smooth function yp in Ai it holds that 

Q{f,^)=[ (—w^-Ky)]F^dsga+f |V(/? --VzpdAgo 
JdM \ OT J Jm ^ 

where r is the outer normal veetor to dAi. Consequently if 

k{y)<^^^, ye dAi 
then A4 is minimizer for the area funetional. 

2 

Proof. The proof follows directly testing equation (12.71) against ^ and integrating by 
parts. □ 
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We directly see from the previous proposition that if Ki , Kq < 0 and Ai is an even 
catenoidal portion with small area, we are in the case nonconvex and M a local mini- 
mizer for the area functional. To see this it is enough to consider the value y > 0 from 
above and Ai to be an even piece of catenoid contained in the portion associated to y where 
z = —Z2 > 0. 

If Ai has large enough area it resembles the entire catenoid which has positive morse index. 
Cutting off an eigenfunction of Am + |^mP associated to a positive eigenvalue in a way that 
the boundary condition is not seen one can obtain a direction where the second variation of 
surface area is negative and therefore a catenoidal portion Ai with large area is unstable. 

The former situation ocurrs also in a general complete embedded minimal surface with 
finite total curvature for which the Morse Index is finite. 


3. Geometrical computations 


In this part we compute the euclidean laplacian in a neighborhood M oi Ai and the normal 
derivate ^ in 511 n AA in system of coordinates that will be suitable for our subsequent 
developments. 

In what follows we consider a large dilation of Ai, denoted b Ma ■= a~^Ai for a > 0 
small. First we compute the euclidean laplacian well inside the set close to Ai. Following 
the developments from mi, for the complete embedded nondegenerate minimal surface M 
with hnite total curvature denote, Mq the part of M inside the cylinder {r(y) < + 1}. 

To parameterize Mq we take a mapping y G C —)• y := yb(y) with associated induced 
metric given by g := {gij) 2 x 2 - laplace beltrami operator of M inside the cylinder can be 
computed as 

Am = (Vdet(y)y*^ (3.1) 

where g~^ = {g^^)) 2 x 2 is the inverse of g and where the coefficient matrix {a^j) 2 x 2 is an smooth 
uniformly elliptic matrix. 


Next, consider the set 

^ := {y = (yi,y 2 ) e : ^(y) > -^o} 


and for M^, the k — th end of M, we use the parametrization 

y G T) Yfc(y) := y^Ci + Fk{Y)e^. 

Next, we notice that the unit normal vector to M at a point y G Mk has the expression in 
coordinates 


^{y) 


(- 1 )" 

Vl + |VFfc(y)|2 


{diFkCi - es) 


(- 1 )^ 63 + afc^ + C>(r 2), y = Yk{y) 


so that diU = 0{r 2) and |AmP = 0{r as r —)• oo. 

In the coordinates Tfc(y) on the metric g := {gij) 2 x 2 satistisfies 
9ij = ^ij + i,j = 1,2 as r oo 

and this relations can be differentiated in the sense that the term gains one negative 

power for every time it is differentiated. 
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We compute the laplace-beltrami operator on using (j.S.lj) to find 

Am = Ay + 0{r~‘^)dij + 0{r~^)di, on M^. (3.2) 

M is parameterized completely by the m + 1 local coordinates described above and we 
observe that expression (j3.1j) holds in the entire M, where for k = 1, ... ,772- the coefficients 
on Mfc satisfy 

bi{y) = 0{r~^) for fe = 1,..., m, as r —)■ 00 . 


We can consider Fermi coordinates given by the mapping X{y,z) := y + zv{y) which 
provides a change of variables in the neighborhood of M 

Af := {x = y + zu{y) : \z\ < y + 6 log (2 + r{y))} 

and where we have the expression 


Ax — dzz + Am — zIAmI"^ dz + D ( 3 - 3 ) 

where the Am is given in ( 13 . 11 ) . ( 13 . 21 ) and where the differential operator D has the form 
D = za\j{y,z)dij + zb]{y,z)di + z^bl{y,z)dz 
and the smooth functions aij{y,z), bi{y,z) satisfy on the ends of M 

141 + = 0{r-% \b}\ + \rVbj\ = 0(r-3) 


(3.4) 


\bl\ + kV 6 ^| = 0 (r-®) 

as r ^ 00 , uniformly on z in the neighborhood M of M, see lemma 2.1 [H]. 

The term 63 comes from the mean curvature of the normally translated surface 

Mz := {y + zu{y) : y £ M} C M 
for fxed 2 : small. It is well known that 


Hmz ■— Hm ~ + ^ 2 ) 4" ^^(^1 4" ^ 2 ) 4" 0 ( 2 ^ r ^*^). 

'-V-' 

where /ci, ^2 are the principal curvatures of M. From the asymptotics of Xm’^, we have that 
ki = 0{r~‘^) and since Hm = ki + k 2 = 0 , it follows that kf + k^ = 0 and therefore the 
expansion for 63 in ()3.4p follows. 

Recall that for small a > 0 Mq, := and let us denote the dilated ends of M by 

Aik,a = ot~^Mk. Next, for a smooth function h defined in A4, we consider dilated and 
translated Fermi coordinates 


Xa,h{y,t) ■= X{y,z), y = ay, z = a{t + h{ay)) 
for y £ Ma and \t + h{ay)\ < ^ + | log (2 + r{ay)). 

Scaling and translating expression (13.3p we obtain 

o? Ax = Ax^t,^ = dtt + Am^ — a^{AMh + \AM\'^h}dt — a^\AM\‘^tdt 

— 2q; a^jdih djt + o? a^jdihdj hdu + Da,h (3.5) 
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where 

Da,h = Oi{t + h)a}j{ay, a{t + h)) {dij - 2adihdit - a^dijhdt + a^dihdjhdtt) 

+ a‘^{t + h)b\{ay, a{t + h)) {di — adihdt) 

+ a^{t + Kfh\{ay,a{t + h)) dt- (3.6) 

Expression (13.51) holds true in the region 

J^a,h '■= a~^9. n iy + (t + h{ay))i'{ay) : \t + h\ < — + - log(2 + r{ay))\ (3.7) 

[ a a ) 

and we will use it to handle equation m well inside the region a ^ (n n Af ). 

In this regard, observe that our geometrical considerations above do not take into account 
the effect of dfl H Af. We need to consider a new system of coordinates to handle boundary 
related computations. 

In this part, we will use that the entire surface M is orthogonal to dffl. Using assumption 
(iii) from the introduction, for fix A: = 1, ... ,m and in we can assume that the closed 
simple curve Ck '■= Mk H U is parameterized by a mapping 

0 G (0,4) i-A 7 fc := 7 fc(u). 

The mapping 7 ^ in turn, allows a smooth orthogonal extension to an open neighborhood 
of Ck in Mk- Abusing the notation we write this extension as 

(p,!;) G^ 7 fc = 7fc(/5,u), pe{-S,6), 6»g(0,4) 

which can be chosen satisfying 

7 fc( 0 ,u) = 7 fc(u), 9„7fe(0,i;) = i9„7fc(u), dpjk T d^'jk 

and such that 7 ([0,(5) x (0,4)) C Mk n U. In essence, the coordinates (p, u) work as polar 
coordinates in Mk near Ck- 

In the coordinates 7 fc(p, u) and ommiting the explicit depedence on k, the laplace beltrami 
operator of M close to Ck, takes the form 

Am = a%dij + b^dj, i,j = p, 6 (3.8) 

and where 

alp{p,0) = |(9p7fc^^ aPpv = ^Ip = 0 - 

f>^p{p^'^) = \dplk\~‘^\dvlk\~‘^ < dpv'yk;d^lk > +|'9p7fcr'‘ < dppjkldp-fk > 

bl{p,v) = \dp'yk\~‘^\dy'yk\~'^ < 5p„7fc;(9p7fc > +|i9p7fcr'‘ < d-^vlk] d^'yk > • 

Proceeding as above and associated to the coordinate system y = 'yk{p,v) we consider 
Fermi coordinates 

Xk{p, V, z) := 7 fc(p, v) + ziy{p, v) 

in the neighborhood of Mk 

^f ■= {'yk{p,v) + zi^ip,v) ■ \z\ <V + 5\og{2 + r{jkip,v))), |p| < 5, 'i;G(0,4)}. 

To described fflCiAf near Ck we assume the existence of a smooth funtion Gk = Gk{v,z) 
such that Gk{v,C) = 0 
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Making a translation of the integral lines of associated to the parameterization jk in 
the p direction using the function Gk, namely 

p{s,v,z) := s + Gk{v,z), |s| < (5, i;e(0,4), \z\ < p 

we can described n AA using the modified Fermi coordinates 

X{s,v,z) := 7 fc(p(s,i;,z),n) + zi'{p(s,v, z),v) 

which actually provide a change of variables in the larger set 

:= jx = X{s,e,z) : jz] < p + 6log{2 + ripk{p{s,v, z), 9))), |s| <6, 9 G (0,4)| 

and clearly 

dQrW = {'yk{Gk{v,z),9)+zi^{Gk{v,z) : 9 G {0,lk), \z\ < p+6log {2 + rijk {Gk{v, z),v)))}. 

Observe that 

d^X{0,v,0) = 1^(0, v) + (9p7fc(0,n) • d^Gkiv,0) 

and 

dyX{0,v,0) =de'yk{0,v), dsX{0,v,0) = dp'yk{0,v) 

so that, by making the inner product of dzX{0,v,0) against (9 sX(0,t;, 0), we observe that 
conditon (I) is rewritten as dzGk{v,0) = 0. 

Summarizing, we assume on the function Gk{v,z) that 

Gfc(n,0)=0, dzGk{v,0)=0, v G i0,lk) (3.9) 

From ()3.9p . the asymptotic expansion in powers of 2 of the mapping X{s, v, z) reads as 

2 3 

X{s,v,z) := 7fc('S,i^) +zv{s,v) + %-qi{s,v) + ^q 2 {s,v) +0{z^) (3.10) 

2 0 

where qi,q 2 -I- u with expressions given by 

qi{s,v) := dp-fk{s,v) • dzzGk{v,0), q 2 {s,v) = dp-fk ■ 5®Gfc(n,0) + 3dpU ■ dzzGk{v,0). 

Taking derivaties in expression (I3.10p and ommiting the dependence on A:, we can compute 
the induced metric of M, which in this coordinates takes the form 

/ |5p7p 0 0\ / -2M -2N <dpr,qi>\ 

~g=[ 0 |5„7|2| 0 +z -2N -2R <dyr,qi> +0(^^) (3.11) 

\ 0 0 1 / \ < dpT, qi> < dyj; qi > 0 / 

where 

—M =< dpj; dpV >, —R =< 9^7; dyV >, —2N =< dp'^] dyV > + < (9„7; dpV > 

where all the entries of the matrices above are evaluated at (s,n). 
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Consequently, the inverse of the metric has the asymptotic expression g ^ = (ff*'^) 3 x 3 has 
the form 

/ |9p7|-^ 0 0 \ 

= 0 |5„7|-2| 0 

Vo 0 1/ 

/ 2 \dp^\-^M ^ 2 \d,^\-^\d^^\fN -\dp^\-^ <d,r,qi> \ 

+ z\ 2 \dp'y\-‘^\dy'y\-^N 2 \dy'y\-^R -|<9„7|"2 < (9^7; > \+ 0 (z^) 

\ -\dpl\~‘^ < dpl^qi > -\dvl\~^ < <9^7;> 0 / 

(3.12) 


We consider again the function h as above and proceeding in the same fashion we take 
dilated and translated Modified Fermi coordinates 

^a,h{s, 9, t) = a~^X{as, a6, a{t + h{as, a9))) 


for 


0 < s < —, 6* G (0, —), \z\ < — + — log(2 + r( 7 (Q;s, a9))). 

a a a a 

After a series of lenghty but necessary computations we arrive to the expression for the 
euclidean laplacian in the coordinates /j 

= Aj^^ ^ = dtt + ^Ma ~ tdt — {Am/i + |^m^} dt + Dq + ^ (3.13) 

where following |39j . we denote 

/i(u) = |(9„7(0 ,u)| > 0, / 2 (u) = |i9p7(0,i;)| > 0, I{v) = hivf dzzGk{v,Q) 


A{v) =< dpp'y{0,v)-,dp'y{0,v) >, C{v) =< dpyj{0,v)-,dy'y{0,v) > 


R{v) =< dyy-/{0, v); (9„7(0, v) >, E{v) =< dpyjiO, u); 5p7(0, v) >, 

to find that 

Do = -2a—-^^^^—{t + h){dst-adphdtt)-al^'^{ad)dyyhdet-al2'^{ae)dpphdst 
ti [o;u )L 2 [cxu j 

Wli^(a9)ldyhl^du + a^l2^(a9)ldphl^du - 2a^sl2^(a9)A(a9)ldphl^dtt (3.14) 


^ ^1 i^^ot A Oi ci2 i^^st G^dphdti) 

+ a%i {dp — adyhdt) + a%2 {ds — adphdt) 

+ a^{t + Kf‘ho{as,a6,a{t + h))dt + a^Ra, (3.15) 

where the functions ai, ai, 6 i, 62 , bs are smooth with bounded derivatives and Ra is a differ¬ 
ential operator having dependence on h and its derivatives. 

We remark that in the coordinates (s, 6) 


^Mr- — 


liiaO) 


dee + 


1 o o -4(a6i) 

ifP) - 2 “*— 


“ 1 “ O! 


E{ae) 


4{a9) 
R{ae) 


+ a 


C{ae) A{ae) 


I\{a6)l2{a9) l^iotO) 


d. 


ll{ae)ll{ae) lf(a9) 


de + O(a^) 


(3.16) 
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Next, we turn our attention to the boundary condition. It can be check directly from (j.S.lip 
and (j3.12ll that the boundary condition reads as 


d 


— = ^/^^p + 


9 


12 


=a. 


9 


13 


--d. 


dn ^ 

so that, after dilating and translating we find that the boundary condition becomes 

ds = al{a9)tdt + a {dph + I{a6)h} dt 


2a{t + h)mi{a9) + a^{t + hYdi{a9) 


ds + 


a{t + h)m2{a9) + a d2{a9, a{t + h)) 
2™ ^a\Pl I u\ I „.3 


dg 


a^{t + hYI{a9)mi{a9) + 2a^mi{a9)dph{t + h) — a'^m2{9)dyh{t + h) + dt (3.17) 
where 

”^ 1 ( 1 ^) = \h{v)\~^ < dpi^{0,v);dp-/{0,v) > 


m2{v) = |/i(u)| ^ (< 5„i/(0,u);(Ip7(0,u) > + < dpv{Q,v)-,dy-f{Q,v) >) 
and Ba = Ba{9, t, h, Vm^) is has dependence in its variables. 

We mention that in the case of the catenoid and an axially symmetric domain, following the 
scheme in [HE] one can parameterize with only one set of coordinates and the calculations 
reduce considerably. 


4. THE APPROXIMATION AND PRELIMINARY DISCUSSION 

The proof of our main result relies on a Lyapunov-Schmitt procedure near an almost solu¬ 
tion to the equation. This section is devoted to find a good global approximation to perform 
this reduction. 

For this, we denote f{u) = u{l — u^) and we consider the solution of the ode 
w"{s) + t(;(l — = 0, s G R, ii;(±oo) = ±1 


which is given explicitly by 


w{s) 



s G R 


and has the asymptotic properties 

w{s) = 1 — -I- O , 

w{s) = -1 + 2e^" + O , 

w'{s) = 2x/2e-^l^l O (^e-2V2|*l) , 


s > 1 
s < —1 
|s| > 1 


(4.1) 


where w' = 

Next, denoting Hq, := and after an obvious rescaling we are lead to consider the 

problem 

du 

Au + f{u) = 0, in fla, —— = 0, on dQ.a (4-2) 

C/TT-ck 

where Ua stands for the inward unit normal vector to dQ,a- 


In what follows for a function U = U{x) and a subdomain 12' C we denote by 


S{U) = AU + f{U), in O'. 
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4.1. The inner approximation. Let us take h £ satisfying the apriori estimate. 

+ II^^IIl°°(A 1) + ll^lll/°°(A^) ^ (4-3) 

where the constant /C is going to be chosen large but independent of a > 0. 

Using the coordinates we set as first local approximation 

uq{x) = w{t), X = Xa,h{y,0,t) G q;“^(0 n AT). 

When computing the error created by uq using (I3.5l) - (|3.6p . in n AA) we find that 

5(tto) = —o^iAMh + \AM\'^h}w'{t) — a'^\AM\‘^t'w' (t) + a^dihdjhw" (t) 

—a^{t + h)a}j{ay,a{t + h)) {dijhw {t) — dihdjhw"(t)) — a^b}{ay,a{t + h))dihw'(t) 

— a'^{t + h)b\{ay, a{t + h))w'{t) (4-4) 

where \AM\,h,dih,dijh are evaluated at ay. 

Observe that if we take h = 0, the size and behavior of the error in expression (14.4p is given 
by 

—a‘^\AMo\‘^tw\t) + a^b\{ay, at)t^w'{t). 

As in [H], due to the presence of the O(a^) term we need to improve this approximation. 
Hence we consider the function 'ipi(t) solving the ode 

dtt'4^i{t) + F 't G IR. (4.5) 

Using variations of parameters formula and the fact that 

f t{w '= 0 

Jr 

we obtain that given by the formula 

rt poo 

=—w{t) / w'{s)~'^ / ^w'd^ds 

Jo Js 

from where it follows at once that 

<Cj, j e , 0<a <V2. 

So, we consider as a second approximation in the region PI A^) the function 

ui{x) = w{t) + (j)i{y,t) (4.6) 

where in the coordinates X^^h 

(f)i{y,t) = a^\AMiay)\‘^ 

When computing the inner error of this new approximation we find that 

S{ui) = A(pi + f{w{t))4>i + S{uo) + f{w{t) + (pi) - f{w{t)) - f{w{t))4>i 

= —a‘^{AMh+\AM\‘^h}'w'{t)+a‘^dihdjhw'\t)—a^{t+h)a}j{ay,a{t+h)) [dijhw'{t) — dihdjhw"(t)) 
—a^bl{ay, a{t + h))dihw'{t) — a!^{t + h)b\{ay, a{t + h))w'{t) 
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+a^AM(|^Mp)i/'i(i) - a'^lAMh + + a^lAMl^dt-ipiit) 

-2a^a^j{ay)dihdj{\AM\‘^)dt'ilJi{t) + a^a^j{ay)dihdjh\AM\^dtt'ipi{t) 

A((/>i) + a^Ri^a{ay, t, h, VMh, Dj^h) (4.7) 

where 

Ni4>i) = f{wit))4>i + S{uo) + f{w{t) + (pi) - fiw{t)) - f'{w{t))(pi ~ 
and the differential operator has dependence of all of its variables with 

|Viii,„| + \Ri,a\ < for 0 < fj < V2. 

From the error (14.7p we see that in PI Af) the open neighborhood of A4 

|5(ni) — + \AMfh}w'{t)\ < (4.8) 


4.2. Boundary correction. It is clear that our approximation ui can be defined in the set 
a“^(fl n A^), but ui does not satisfy in general the boundary condition. In this regard we 
need to make a further improvement of the approximation ui by adding boundary correction 
terms. 


Let us consider a cut-off function /3 such that 


Pip) = 


1, 0 < p < I 

0, p > 6 


For the k — th end of a ^M, we consider a cut-off function 13a = Pk,aix) = /3(as) for 
X = Xa,his,0,t). 

Near the boundary we consider a new approximation of the form 

m 

U 2 ix) = Ui{x) -F ^ iPa,k4>2,ki^) + Pa,k4>3,ki^)) 
k=l 

where (p 2 ,k, 4’3,kix) are going to be chosen of order 0{a) and 0(0^) respectively. To see how 
to choose 4 > 2 ,k, (p 3 ,ki we proceed first by computing the error of the boundary condition created 
by the approximation U 2 using the coordinates Xa^h and expression (|3.17l) for a fixed end ■ 
In what follows we ommit the explicit dependence of k but we remark that the developments 
in this part hold true regardless of the end we are working with, since they all have the same 
geometrical structure and the supports of the cutt-off functions 13a,k within te region Afa,h 
close to every end Mk,a are far away from each other. 

Recalling that we have assumed h = Oy/ 2 ,p(^j^'^[a) and splitting the boundary error in 
powers of a we find 

B{u 2 ) = -ds(l )2 + al{a6)tw'{t) - ds4>3 


+CX {dphi -|- I{ad)hi) w'{t) — a^I{a9)mi{a6)t^w'{t) -|- al{a9)tdt(p2 — 2ami{a9)tds(p2 
-a^dp{\AM\^)ipi{t) + a^I{a9)\AM\'^dt'4)i{t) 
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+a {dph + I{a9)K} dt<p2 — 2mi{a6)hds4>2 — o?{t + di{a6)ds(l>2 

— 2a^ I{a9)mi{a0)hw' {t) — I{a9)mi{a9)t^dt(j)2 + (4.9) 

where the term -Bo,a satisfies that 

\VBo,a\ + \Bo,a\ 

Our goal is to get a boundary error of order Cl(Q;^e“'^|t|). In order to do so, we hrst choose 
4)2 solving the equation 

+ Am„ 02 +/'(w^(t))02 = 0, in Mq, X IR 
ds42 = al{a9)tw'{t) 

from where we obtain that is odd in t and from proposition 16.11 it follows that in the 

norms dSSIl-dOl) 

\\D‘^4>2\\p,a + \\e‘^^^^'^4>2\\L°°{McXR) + l|e'^'*'02||L°°(Mc<XR) < Ca. 

To choose (/> 3 , we make the decomposition 

t‘^w'{t) = ciw'{t) + gi{t), Cl = \\w'\\l2(^^) j t^{w'{t)fdt 

J R 

tds(t)2i{-,t) = C2{-)w'{t) + g2{-.,t), f g2i-,t)w'{t)dt = 0 

Jr 

and we write the second line in (14.9p as 

+a {dph + I{a9)h) w'{t) — a^cil{a9)mi{a9)w'{t) — 2a'mi{a9)c2{9)w'{t) 


—a^I{a9)mi{a9)gi{t) — 2ami{a9)g2{9, t) + al{a9)tdt(f)2- 

Hence to choose (ps we hrst ask (p 3 to satisfy the boundary condition on dM^ x R 

ds4>3{0,9,t) = —a^I{a9)mi{a9)gi{t) — 2ami{a9)g2{9,t) + al{a9)tdt42- (4.10) 

Next we compute the error of the approximation U 2 near the boundary using expression 
()3.13p . Setting 

Sout{u2) = S{ui) + A02 + f{w{t))42 + A(/>3 + f'{w{t))43 + ^f'\w{t)){42 + (psf 
[fiui) - f{w{t))]{42 + h) + - f{w{t))]{42 + 

[f{ui +42 + 43 ) - f{ui) - f{ui){42 + 43 ) - ^f''{ui){42 + 43 )^]- 


So that we can write this error explicitely in coordinates (s, 9, t) as 

Sout{u2) = S{ui)+dtt43+^Mc,43+f'{w{i))43-oP‘\AM\'^tdt42+\f''{w{t))4l-2a nJ dst42 

2 Ii[a9)l2{a9) 

-c? {AmH + dt42 - 2a—^^^p—^hdst42 + ‘2-c? , a\ + h)dphdtt42 

Li[a9)L2[a9) Li[a9jL2{a9) 
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-al^ ‘^{a9)dvvh^dgt4>2 - al2 ‘^{a9)dpphdst4>2 + oi^li {a9)\d^h\^dtt4>2 
+a‘^l2‘^{a9)\dpph\‘^dtt4>2 - 2a^sl2'^{a9)A{a9)\dph\‘^dtt4>2 
+a‘^di{as, a9, a{t + h)) {det4>2 - adyh(j)tt4>2} + a^ 02 (as, a9, a{t + h)) {dst(l >2 - oidph4>u4>2} 


+a‘^bi{as,a9,a{t + h)) {dg(f )2 — advh(j)t(j) 2 } + a‘^b2{as,a9,a{t + h)) {ds4>2 — ocdph4>t4>2} 

+f'{w{t)){(/)2-(j)3 + ^(j4) + [f'{ui)-f'{w{t))]{(j)2+(l)3) + +a^{t+h)b3{as,a9,a{t+h))dt(j)2+R2,a 

(4.11) 

where i?2,a = R2,a{ots^a9^t,h,'S/j^h^D'j^h) = 0{a‘^) and 

\DR2,a\ + 1-^2,al < Ca^. 

we notice that 

5(1X2) = (1 - /3a)5(ltl) + PaSoutim) + 2Vl3a ' 4>2 + 2VjSa ' V^s + {4>2 + 03)A/3a 

/(ill + Paih + 4>3)) - f{ui) - Pafiui + h + (f^s) “ (1 “ /3a)/(«l) 
and observe that 

V/3„ • V02 = adsp{as)ds(t)2 = (4.12) 

and 

/(ill + I3a{4>2 + (t>3)) - f{ui) - Pafiui + 4>2 + (p3) “ (1 “ Pa)f{ui) = 

Hence in order to improve the approximation we need to get rid of the terms in the hrst 
line of (I4.1ip and the term in (j4.12p so that we choose 03 solving the linear problem 

dtt4>3+^Mc,4>3+f{w{t))4>3 = a'^\AM\‘^tdt(t)2-l:f'{w{t))(l)l+2a—^^^j-^tdst4>2-ads(3{as)ds(t)2, in Mq,xR 

2 Li[a9)L2[a9) 

ds(p3{0,9,t) = —a^I{a9)mi{a9)gi{t) — 2ami{a9)g2{9,t) + al{a9)tdt4>2- 

/ 03(-, t)it>'(t)(it = 0, in Mq. 

Jr 

From this, 03 satisfies 

||-D^03||p,o- + V0311XR) + l|e'^'*'</’3||L°°(Mc<xR) < Ca^. 

From expression ()4.1ip we directly check that 

5(112) = 5(iii) + Eo,„ + ^i,„ (4.13) 

where Ei^a = Ei^a{as,a9,t, h,VMh, Dj^h) = 0(a^+*) and 

\RRi,a\ + ^ 1 = 0 , 1 . 

and from (14.91) and (I4.10p the boundary error takes the form 

B{u2) = OL {dph\ + I{a9)hi) w'{t) — cil{a9)Tni{a9)w'(t) — 2ami(a0)c2(0)ii''(t) (4.14) 


Bi^a + H 2 ,, 
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where \\c2\\L°°{dMa) ^ C'a and = Ei^a.{as, aO, t, h, Vj^h, Dj^h) = 0{a^ + i) and 

\DBi^ci\ + i = 1,2. 

To get the right size of the boundary error we impose on h the boundary condition 
dph + I{v)h = acil{v)mi{v) + 2 mi(u)c 2 
where we remark that ||c 2 < Ca so that 

\\acilmi + 2miC2 \\L°°{dM) ^ Ea. 


4.3. The global approximation. Observe that so far, our approximation is defined only in 
the open set f^N'). In order to define a global approximation in Oq, let us define the 

function 


f +1, x G a 
\ —1, X G a~^S~. 


(4.15) 


which is clearly an exact solution of the equation wherever it is smooth. 


The idea to get a global approximation is to consider the approximation U 2 well inside 
while outside Ma,h we interpolate with the function (x). In order to make this precise, 
let us take a non-negative function f3 in C'°°([R) such that 

Pis) = { 1, |s| < 10, |s| > 2 

and consider the following cut-off function in N'a,h given by 

/3^(x) = Pi\t + h{ay)\ - - - <5ln(r(a?/)) + 2), x = Xa,h(y,^,t) S Afa,h- 

a 


With the aid of this, we set up as approximation in Qa the function 
[/(x) = Pp(x)u 2 (x) + (1 - Pri(x) (x), X G fla- 
and we compute the new error created by this approximation as follows 


(4.16) 


S(17) = AU + f{U) = P2 S{u2) + E 

where 

E = fiPnU) - pr,fiU) + 2VPp • Vu + uA/3^. 

Using that z = \t + h(ay)| we see that the derivatives of /3i do not depend on the derivatives 
of h. On the other hand, due to the choice of Pr] and the explicit form of E, the error created 
only takes into account the values of P^ in the set 

X = X^ hiy,G,t) GJ\fa,h, \t + h{ay)\ > - + 4:ln{r{ay)) - 2, 

a 

so we get the following estimate for the error E 

\E\ < Ce-i. 
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5. The proof of Theorem [T] 

The proof of Theorem [T] is fairly technical. So, to keep the presentation as clear as possible, 
we sketch the steps of the proof and in the next sections we give the detailed proofs of the 
lemmas and propositions mentioned here. 

First, we introduce the norms we consider to set up an appropriate functional analytic 
scheme for the proof of Theorem [H For q;> 0, l<p<oo and a function /(x), defined in 
Qa, we set 


WfWp,^ := sup \\f\\LP{Br{x))- (5-1) 

We also consider for functions g = g{y,t), (j) = defined in the whole Mq, x R, the 

norms 

\\9\\p,a ■= sup (5-2) 

{y,t)eMaXR 

||</>||2,p,<t := + \\D4>\\oo,a + ||</>||oo,o-- (5.3) 

where dA^ ■= dyg^ dt. Of course, in the case p = +oo, we have that L°°{Bi{y, t)), dA^) = 
L-(Si(y,t)). 

While for a function G defined in dMa x R we consider the norm 

l|G'llp,fT := l|e‘^'*'G'||x,p(aMcXR) 

And we mentioned in the previous section, we set the norm for the parameter h as 


ll^ll* — \\D'^h\\LP{M) + II^^IIl°°(AI) + II^IIl°°(A^) (5-4) 


We look for a solution to equation (|l.lll of the form 

Ua{x) = U{x) + (f{x) 

where U{x) is the global approximation defined in (I4.16P and ip is going to be chosen small 
in some appropriate sense. Thus, we need to solve the problem 

Ap + f{U)p + S{U) + N{p) = 0 


or equivalently 


+ fiU)p 


-S{W)-N{p) 
-l3pS{u2) -E- N{p) 


(5.5) 


where 


N{p) = fiU + p)- f{U)- f'{U)p. 
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5.1. The Gluing Procedure. In order to solve problem (15.5h . we consider again the cut-off 
function /3, from previous section, and we define for every n G , the cut-off function 


Cn{x) : — 


_/ I3{\t + h{ay)\ - ^ + n) if x = Xa,hiyi,y 2 ,t) E J\fa,h 


0 


if X e Ma,h- 


(5.6) 


We look for a solution to (j5.5ll with the particular form 

y^{x) = C,2{x)(t){y,t) + '4){x) 

where (j){y,t) is defined for every {y,t) G M„ x IR and tplx) is defined in the whole flo,. So, we 
find from equation (15.511 that 

C2 + f{U)cl> + CiU^ + S{U) + C2iV(0 +1/.)] 

+ AiA - [2 - (1 - C2)[f'{U) + 2]]^|; + (1 - C2)SiU) 

+ 2 VC 2 • ^ Na,h4> + 0A^2 + (1 ~ C2)a[^20 -i- V'] = 0. 

Hence, we will have constructed a solution to the problem (|5.5I1 . if solve the system 

+ C2f^'0 + ^(C/)-RC 2 A((/>-I-V') = 0, in |t-R/i(a 2 /)| < - - 1 (5.7) 


a 


AV^ - [2 - (1 - C2)[nu) + 2]]iA + (1 - C2)S{U) 

+ 2 VC 2 • + <(>AC 2 + (1 - C 2 ) A[C20 -t- lA] = 0, in He 

As for the boundary conditions we compute 

dC,2 , '9V’ 


(5.8) 




du2 , ^ d(f) ^ ^ 


dua duoi dua ^ dUa ' OUa 


+ 


= 0 . 


Therefore, as we proceeded above, we reduce the boundary condition to the boundary 
condition system 


du2 d(t> 

+ C2t^— - 0 


dUr 


‘drin 


5V’ , e ^ , ..^^2 


= 0 . 


(5.9) 

(5.10) 


Next, we extend (I5.7jl to a qualitative similar equation in Mq. x R. Let us set 

Ri(f)) := C4[Aa4,,^ - dtt - Am„]- 

Observe that R{4>) is understood to be zero for \t + h{ay)\ > ^ + 2 and so we consider the 
equation 

du4> + Amc<V + = -S{u2) - R{4>) 


- {f'{u 2 ) - f{w{t)))(j) - C2tt2V’ - C 2 N{(f) + 'ljj), in Ma X R. (5.11) 

where from expression (|4.1ip and ommiting the depedence on k, we have on every the k — th 
end Mfc Q, explicitely 

S{u2) = S{ui)-a'^ {AmH + |AMph} (9tV’2-2« /n9giV*2+2«^ ^J f+/?-)5ph5ttV*2 

Li[aa)L2[aU) li{at>)l2[ad) 

+2adp/3{as)ds(j)3+Oi^dpp/3{as){(p2+4>3)-ctli^ ia6)d^,^,h'^ dgt(i>2-al2‘^ ict0)dpphdst(p2+a^li‘^ iaG)\dvh\'^ dtt4)2 
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+a^l2'^{ae)\dpph\‘^dtt(j)2-‘^a^sl2'^iae)A{a9)\dph\‘^dtt(i>2+f'{w{t)){4i2-(l}3+^4>l)+[f{ui)-f'{w{t))]{4i2+(p3) 
+a‘^C 4 :aiias,a 6 ,a{t + h)) {dgt 4>2 - adyh(j)tt 4 > 2 } +( 4 ^ 2 ( 0 ( 3 , 09 , a{t + h)) {dst 4>2 - 0 (dph(t)tt(l> 2 } 


+a^(4bi{as,a9,a{t + h)) {dg4>2 - advh4>t4>2} + CJ> 2 {c(S, a9, a{t + h)) {ds(f )2 - adph(j)t4>2} 


+ a^C 4 (t + h)h'i{as, a9, a(t + h))dt4>2 + ( 4 ^ 2,0 (5-12) 

and from expression (14.71) we write 

S{ui) = —a^{AMh + \AM\'^h}w'{t) + a^dihdjhw" (t) 


AM{\AM\‘^)'tpi{t) — a^{AMh + + ot^\AM\^tdt'tpi{t) 

-2a^a^i^(ay)dihdj(\AM\^)dt'4^i(t) + a^a^ij(ay)dihdjh\AM\‘^dttil2i(t) 
—a^C,i(t + h)a\j(ay,a(t + h)) (dijhw'(t) — dihdjhw"(t)'j 
—a^C, 4 b\(ay,a(t + h))dihw' (t) — a^(t + h)C,ih\(ay,a(t + h))w'(t) 


N((j)i) + a^URi^aioy, t, h, VmK D\jh) 


(5.13) 


Observe that S(ui) and 5( 1 x 2 ) coincide with S(ui), S'( 1 x 2 ) but the parts that are not dehned 
for alH G IR are cut-off outside the support of ^ 2 - 

We proceed in the same fashion for the boundary condition we and writing 

d 


13 = (4 


V9^ 


dUr 


- ds 


it suffices to consider cj) satisfying 

dr^(l) + B{(j)) = B{u2) 

where Tq = s is the tangent inward direction to dMa and in expression (14.91) we cut-off the 
parts that are not dehned for every t. We write also for further purposes 

B{u 2) = -a^dp{\AM\^)il’i{t) + a^I{a9)\AM\^dt'il2i{t) 


+a {dph -|- I(a9)K} dt4>2 — 2mi(a9)hds(l)2 — 0^(1 -|- Kf (4,d\{a9)ds4>2 

— 2 a^ I (a9)mi(a9)hw' (t) — I (a9)mi(a9)t‘^ dt(t >2 + (5-14) 

Observe that again we have ommitted the depedence on the end Mk^a for notational con¬ 
venience. 

We solve hrst (I5.8l) - (|5.10p . using the fact that the potential 2— (1 —C2)[/'(C7)+2] is uniformly 
positive, so that the linear operator behaves like A — 2. A solution x/^ = is then found 
from the contraction mapping principle. We collect this discussion in the following lemma, 
that will be proven in detail in section 5. 
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Proposition 5.1. Assume and 3 < p < oo and let h be as in (I4.3p . Then, for every a > 0 
sufficiently small and every cf such that ||<^|| 2 ,^, 0 - < 1? equation dSS]) has a unique solution 
if = Even more the operator '!'(</>) turns out to be lipschitz in cj). More precisely, 

satisfies that 

HWx ■= + ll-C>i/>||oo + llV’lloo < (5.15) 

and 

CTl 

||4'((/)i) - ^{4>2 )\\x < Ce~~\\(l)i - (l) 2 h,p,a- (5.16) 

Hence, using Proposition 15.11 we solve equation (|5.1ip with = 4'((/>). Let us set 


N(0) 

:= Ri4>) + ifim) - f'{w{t))(p + C 2(^^2 - (a;))^(</>) + C 2 N{ 4 > + T((/>)), 

in Ma X IR. 

So, 

we only need to solve 



dtt4> + + f{w{t))(l) = - S{u 2 ) - N((/)) + c{y)wfit), in x IR. 

(5.17) 


= B{u 2 ), on dMa x IR 

(5.18) 


/ (f){-,f)w'{t)dt = Q, y ^ Ma 

Jr 

(5.19) 


To solve problem (|5.17p - (l5.18p - (l5.19p . we solve a nonlinear problem in cj), that basically 
eliminates the parts of the error, that do not contribute to the projections. 

The linear theory we develop to solve problem (I5.17p - (l5.18p - ()5.19p . considers right hand 
sides and boundary data with a behavior similar to the that of the error S{u 2 ) and B{u 2 ), 
that as we have seen, is basically of the form 0(e“'^l*l). 

Using the fact that N((^) is Lipschitz with small Lipschitz constant and contraction mapping 
principle in a ball of radius O(a^) in the norm || • II 2 ,p,o-) we solve equation ()5.17p - ()5.18p - ()5.19p . 
This solution 4>, defines a Lipschitz operator (p = 4>(h). This information is collected in the 
following proposition 

Proposition 5.2. Assume 3 < p < 00 and a > 0 is small enough. There exists an universal 
constant C > 0, such that (|5.17l) - (l5.18p - (|5.19l) has a unique solution 0 = <b(/i), satisfying 

Uhp,. < Ca^ 

and 

\mhi)-Hh2)\\2,p,a<Ca^h-h2\U. 


5.2. Adjusting h, to make the projection equal zero. In this part we set cq 
T o conclude the proof of Theorem [H we adjust h so that 



<S'('«2) + N((/)) w'{t)dt = 0. 


Let us integrate (I5.17p - (l5.18p - (j5.19l) against w'{t) to find that 


w 


||2 


S{u 2 )w'{t)dt = {l-fia) [ S{ui)w'{t)dt+ f S{u 2 )w'{t)dt+OL<^(^M^){a‘^) 

Jr Jr 


B 


clearly fia does not depend on t so that we can compute from (I5.13|) 
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—a^ / C)4(i + h)ajj{ay,a{a{t + h))) \^dijhw'{t) — dihdjhw"(t)'^ w'{t)dt 

Jr 

—a^ / C 4 {t + h)b\{ay,a{t + h))dih{w '/ {t + h)^( 4 bl{ay,a{t + h)){w'{t))'^dt 

Jr Jr 

+0^174^1^ / tdt'ipi{t)w'{t)dt + a^Pi{ay,h,Vh,D'^h). 

Jr 

Next we write the reduced error near the boundary. From (j5.12p we obtain that 

[ S{U2) = [ S{ui) - f tdt4)2w'{t)dt-2a-—^^^^p-—h [ dst(j)2w'{t)dt 

Jr Jr Jr hiae^hiaO) 


+2a 


R Jr Jr Il{a0)l2(a0) jp, 

2 


Ii{a9)l2{a9) Jr 


dp / {t+h)dtt4>2w\t)dt—al^‘^{a9)dyyh / d0t4>2'w'{t)dt—al2'^{a9)dpph / dst4>2w'{t)dt 

Jr Jr Jr 


+a^ / C,Ab-i{ots,a9,t){dQtJ)2-otdydtt<p2}w'{t)dt+a^ / C,4d2{as,a9,t) {dst<i)2 - OLdpdtt(t>2} w'{t)dt 

Jr Jr 

-\-a^ / CJ)i{as^a9^t){dQ(j)2 —OLd^dt(l)2}‘^\t)dt+a^ / C^J)2{as^a9^t){ds(l>2 —OLdpdt(l)2}'^\^)dt 
Jr Jr 


f C4{t+h)bl{as,a9,a{t+h))dt(j)2w'{t)dt-2a^—-^^^^—{t+h) f {dst(l>3 - adppduh} w'{t)dt 
Jr h[a9)l2{a9) Jr 


[f'{ui) — f'{w{t))]{(j )2 + 4 >3)w'{t)dt + a^i?o,a(as, ct9) + a‘^ Ra{oiS, a9, h, V/i, D‘^h)w'{t)dt. 


Also observe that from the nonlocal terms and condition (j4.3jl we have 

q{ay,h,Vh,D‘^h) = ['N{(j))w'{t)dt, \\Qi-, h,Vh,D‘^h)\\LM > Ca'^~p 

Jr 

Which implies that 

a~‘^ J (^S{u 2 ) + N((/))^ w'{t)dt = —co{AAf/i+|AAf|^/i}+aPo(y) h,Vh, D‘^h)+a^~pPi{y, h,Vh, D'^h) 

(5.20) 

where 

||A|| + ||/’o|| <C^. 

\DPo\ + \DPi\<C. 

As for the boundary condition, directly from (j5.14p we ask h to satisfy 

dph + I(a9)h = acil {v)mi{v) + 2mi(r;)c2(—) (5-21) 

a 

where we recall that 

\\c2\\L°°(dMc) < Ca 

and the right hand side in (j5.2ip does not depend on h. 

We solve then 

co{Am/i+ |AmP/i} = aPo{y,h,Vh,D‘^h) + a^~pPi{y,h,Vh, D‘^h) 


(5.22) 
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with the boundary condition (j5.21ji as a direct consequence of the theory developed in section 
[2] and a fixed point argument for h in a ball or order 0 (a) in the topology induced by the 
norm|| • ||*. This completes the proof of our theorem. 


6 . PROJECTED LINEAR PROBLEM 


In this part we provide the linear theory for the problem 

dtt(t> + Amc</> + f{w{t))(j) = g + c{y)w'{t), in Mq, x R 

(bw'dt = My), y G ^ = h, on dMa x R 

OTr, 


L 


( 6 . 1 ) 

( 6 . 2 ) 


relies strongly on the fact that solutions to 

dtt(l) ++f{w{t))(t) = {), in X R 

^ = 0, on dMa X R. 

QTq, 

are the scalar multiples of w'{t). The proof folios the same lines of lemma 5.1 in [TT]. We 
simply remark that when decomposing the solution (p as 

(p = c{y)w'{t) + M 

from maximum principle one obtains that \M{y, t)| < for somoe 0 < u < \/2- Defining 

My)= [ 

Jr 

it follows that for certain positive constant A 


-Am^iP + XpJ <0, 


dip 

dTa 


= 0 


where Ta is the iward unit tangent to dMa in M^. Clearly it follows that ip = Q and conse¬ 
quently c{y) is a constant function. 

Proceeding as in section 3 in [9] it suffices to solve the case A = /i = 0 and JgW'{t)dt = 0. 
To prove existence we set 

< p,'tp >:= / Vp-Vip+ 2(p ■ tp 

Jm^xr 

and we consider the space H of function p G H^{Ma x R) such that 

p ■ w' = 0. 


I 

Jh 


MtvXR 


Since f'{w{t)) = —2 -I- 0(e ^1*1) as |t| ^ oo, the equation can be put into the setting 

{I + K)P = f, inH 

where K : H ^ H is a compact operator. From fredholm alternative it follows that existence 
and 

II<(’IIl 2(M„XR) < C'llfl'llL2(M„XR)- 

As for the apriori estimates, we can proceed using a blow up argument following the same 
lines as in the local elliptic regularity as in m- We remark that in our case we also need 
to consider two limiting blow up situations the case of R^ x R when taking limit well inside 
Ma X R and the case of the half space R+ x R^ when taking the limit in coordinates close to 
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dMa X R. The former case is reduced to the case of x IR as limiting situation by using an 
odd reflection respect to the boundary of R+ x R^. 

Thus we have proven the following proposition 

Proposition 6.1. For every p < 3 and for every a > 0 small enough and given arbitrary 
functions g defined in x R and G defined in dMa x R such that 

||S'||p,(T + ||G||p^(j < oo 

there exists a unique bounded solution 4> to problem (I6T])-(I621) satisfying the apriori estimate 

+ ||T>(/)||oo,(T + ||0||oo,o- < C'(llfi'llp,cr + IIG'llp^cr) 

where the constant C dependens only on p > 0. 

7. GLUING REDUCTION AND SOLUTION TO THE PROJECTED PROBLEM. 

In this section, we prove Lemma 15.11 and then we solve the nonlocal projected problem 
(I5.17h - (l5.18l) - (j5.19p . The notations we use in this section have been set up in sections 0] and 

El 


7.1. Solving the Gluing System. Given a fixed (f such that ||(/!)|| 2 ,p ,(7 < 1; we solve problem 
()5.8I) with boundary condition (|5.10l) . To begin with, we observe that there exist constants 
a < b, independent of a, such that 

0 < a < Qaix) < b, for every x G R^ 


where Qa{x) = 2 — (1 — C 2 )[f'{U) + 2]. Using this remark we study the problem 

Aif - Qaix)^^ = g{x), xeita 


= G{x), on 

for given g, G. Concerning solvability of this linear problem we have the following lemma. 


(7.1) 


Lemma 7.1. Assume 3 < p < oo . There exists a constant G > 0 and oq > 0 small enough 
such that for 0 < a < oq <^'^d any given g, G with 

llfl'llLP(Oa) + l|G'llLP(aD<,) < oo 

equation has a unique solution ^|J = ifig), satisfying the a-priori estimate 

Mix < C'(||fl'llLp(ou + IIG'llLp(aUc)) 

The proof of this lemma is standard and we refer the reader to section 2 in [9] for details. 
Now we prove Proposition 15.11 Denote by X, the space of functions if G lU^’^(Dc) such 
that IIV’IIa < oo and let us denote by r(g', G) = ip the solution to the equation (I7.ip . from the 
previuos lemma. We see that the linear map T is continuous i.e 

||r(<7,G')||x<G(||5||LP(o„) + ||G||iP(ao.)) 

Using this we can recast (|5.8p as a fixed point problem, in the following manner 

V' = -r ((1 - C 2 )S{U) + (1 - C 2 )N[C 2 cP + iP], iu 2 - (x))|^ + (7.2) 

Let us take cp, /isatisfying that 


||<(>||2,p,f7 < 1, ll^ll* < Ka. 
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Estimating the size of the right-hand side in (I7.2p . 

Recall that S{U) = Q 2 S{u 2 ) + E. So, we estimate directly to get 

1(1 - C2)S(n2)| < 1*1(1 - C 2 ) < ^ 

this means that 

|(l-C2)5(n2)| 

and so 11(1 - C 2 )S(C/)|Up(o„) < Ca\-^l. 

As for the second term in the right-hand side of (17.21) . the following holds true 

|2VC2-V(/. + 0AC2| < C(l-C2)e-'^l'l||</>||2,p. 

This implies that 

||2VC2-V</. + </>AC2||oo<Ce-'=^. 


Proceeding in the same fashion, we clearly obtain that the boundary condition satisfies 


, ±9C,2 


{U2- {x))^ + 


< Ce 


dua dUc 

Finally we must check the lipschitz character of (1 — C 2 )A*[C 2 </' + i/’]- Take V'lj "02 £ A. Then 

1(1 -C2)A*[C2(/' + V’i] - (1 -C2)A*[C2(/' + 'i/’2]| < 


< 

(l-C2)|/(t/ + 

020 + 0l) 




+ 

+ 

1 

02)-/'(C/)(01- 

-02)1 


< 

C7e-'^^(1 - C 2 ) 

sup |Cl0 + i01 

+ (1- 

0021101 



ie[o,i] 



< 

ce-^^m\oo,. 

+ ||0i||oo + II02I 

|oo)101 

- 02| 


So, we see that 

11(1 - C2)A*[C20 + lAl] - (1 - C2)A*[C2(/> + '0i]||oo < IIV'l - -021100 

In particular we see that ||(1 — C2)A*(C20)||oo < C'e“‘^«. Consider T : A —)• A, T = r(i/^) the 
operator given by the right-hand side of dZS]). From the previous remarks we have that T is a 
contraction provided a is small enough and so we have found ^|J = r(0;) the solution to (I5.8p . 


We can check directly that 'I'(0) = -0 is lipschitz in 0 i.e 


||T(0i) - T(02 )||x < C||(l - C2)[iV(Cl01 + 'k(0l)) - A(Ci02 + 'k(02))]||oo,M 

-|-C'“°'~||01 - 02||2,p,o- 

< Ce-'^i (||T(0i) - T(02)|U + 1101 - 02||2,p,<.) 

Hence for a small, we conclude 

||T(0i)-T(02)||x <C'e-^“||0i -02||2,p,<x. 
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7.2. Solving the Projected Problem. Now we solve problem (I5.17h - ()5.18p - (l5.19l) using 
the linear theory developed in section (U together with a fixed point argument. Prom the 
discussion in 17.11 we have a nonlocal operator ijj = 

Recall that 

N{(p) := R{(p) + (/'(u2) - f'{w{t))(l) + C,2{u2 - (a^))^( 0 ) + C,2N{(f) + {(j))) , in x IR. 

Let us denote 

Ni{(t)) := R{4)) + [f{u2) - /'('w(t))] 0 

N2{4>) ■■= ( 2(112 - {x))^{4>) 


iV3(</>) :=C2N{(j) + ^{(j))) 

We need to investigate the Lipschitz character of W, f = 1, 2,3. We see that 

- N3{^2)\ = C2|iV(-/>l + ^(0l)) - N{(^2 + ^(02))| 

< CC 2 sup + (1 - T){(j )2 + 'I'(02))| • 101 - 02 + ^(0l) “ ^(02)| 

Te[0,i] 

< C [|^'(02)| + 101 - 02| + I^(0l) - ^(02)1 + I02|] ' [|01 “ 02 | + |^(0l) “ ^(02)|] • 

This implies that 

I|-^3(0l) — -^3(02)||p,cr < C[e + ||0l||p,cr + ||02 |Ip,(t] ' ||01 “ 02||p,o-- 

Now we check on Ni(0). Clearly, we just have to pay attention to 7?(0). But notice that 
ii(0) is linear on 0 and 


-2ah'{ay)dty(j) + a'^[h'{ay)fdtt4> + Da^h{(t>)- 
Hence, from the assumptions made on h, we have that 


dt4> 


||iVi(0i) - iVi(02)||p,cr < C'a||0l - 02||2 ,p,(7- 


Observe also that under the assumption made on h we have 

||*5(ii2) + Q:^{AMh + \AM\'^h}w (t)||p,(T < Ca^ 

Hence we for ||0||2,p,o- < Aa^ we have that ||N(0)||p^o- < Ca^. 

As for the boundary condition we check directly from expressions (|3.17p . (I5.18P and ()5.18p 
thta on every end „ the following estimates hold 

\\B{u2)\\oo,a < Ca^ ||H(0)|U,. < Ca(||V0||oo,. + I|0||oo,a) 

with H(0) linear in 0. 

Setting T{g,G) = 0 the linear operator given from proposition 16.11 we recast problem 
()5.17p - ()5.18l) - (j5.19p as the fixed point problem 

0 = T{-S{u 2 ) - N(0), Biu 2 ) - H(0)) =: r(0) 

in the ball 

:={0eX/||0||2,p,.<Aa3} 

where X is the space of function 0 G x IR) with the norm ||0||2,p,(7- Observe that 

||r(0i) - r(02)||x < C'||N(0i) - N(02)||p,. < Ca||0i - 02||x, 0 e B^. 
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On the other hand, because C and A are universal constants and taking A large enough, we 
have that 

\\T{ct>)\\x < C(||5(u 2)IU + l|N(0)IU + \\B{u 2 )\\oo,. + mmoo,.) < Aa^ ct> G . 

Hence, the mapping T is a contraction from the ball B^ onto itself. From the contraction 
mapping principle we get a unique solution ip as require. We denote the solution to (I5.17P - 
()5.18l) - (|5.19l) for h fixed. 

As for the Lipschitz character of <h(/i) it comes from a lengthy by direct computation. We 
left to the reader to check on the details of the proof of the following estimate 

mhi)-^{h2)h,p,a<Ca^\\hl-h2\U. 
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